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Another thread in the tapestry of stellar feedback: X-ray 
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ABSTRACT 

We consider X-ray binaries (XBs) as potential sources of stellar feedback. XBs ob- 
servationally appear able to deposit a high fraction of their power output into their 
local interstellar medium, which may make them a non-negligible source of energy 
input. The formation rate of the most luminous XBs rises with decreasing metallicity, 
which should increase their significance during galaxy formation in the early universe. 
We also argue that stochastic effects are important to XB feedback (XBF) and may 
dominate the systematic changes due to metallicity in many cases. Large stochastic 
variation in the magnitude of XBF at low absolute star formation rates provides a 
natural reason for diversity in the evolution of dwarf galaxies which were initially al- 
most identical, with several percent of such halos experiencing energy input from XBs 
roughly two orders of magnitude above the most likely value. These probability distri- 
butions suggest that the effect of XBF is most commonly significant for total stellar 
masses between approximately 10^ and 10^ M©, which might resolve a current problem 
with modelling populations of such galaxies. We explain how XBs might inject energy 
before luminous supernovae (SNe) contribute significantly to feedback and how XBs 
can assist in keeping gas hot long after the last core-collapse SN has exploded. Energy 
input from XBs produces different behaviour to that from SNe, partly since the peak 
energy input from a mean XB population continues for « 100 Myr after the start of 
a starburst. XBF could be especially important to some dwarf galaxies, potentially 
heating gas without expelling it; the properties of XBF also match those previously 
derived as allowing episodic star formation. We also argue that the efficiency of SN 
feedback (SNF) might be reduced when XBF has had the opportunity to act first. In 
addition, we note that the effect of SNF is unlikely to be scale-free; galaxies smaller 
than « 100 pc might well experience less effective SNF. 
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1 INTRODUCTION 

A collection of processes jointly known as "feedback" seem 
to be vital in governing the state of the baryonic content of 
dark matter haloes; understanding these processes is there- 
fore extremely important to the study of galaxy evolution 
and structure formation. The current-day baryon content 
of an L* galaxy's halo seems to be lower than the cos- 
mic mean, suggesting that baryons have been ejected from 
them, and feedback appears to have an even greater effect 
on the lowest-mass and highest-mass galaxies, truncating 
the high-mass end of the galaxy mass function and flat- 
tening the low-mass end (see, e.g., iKaufFmann et al.lll999l : 
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Benson et al.l l2003l: ICole et all 1200 ll: iTegmark etHI 12004 
Croton et al.ll2006l : iBower et al.ll2006l ). 



Although it is agreed that feedback mechanisms are 
required to produce the properties of galaxies in the ob- 
served Universe, the details of feedback remain far from re- 
solved. It is widely accepted that stellar feedback is dom- 
inant in low-mass galaxies, the canonical sources of which 
are supernovae (SN) and irradiation and winds from massive 
stars. The long-known uncertainty in SN feedback (SNF) 
has been the efficiency with which supernova ex plosions 
can c ouple to the interstellar medium (see, e.g.. iLarsonl 
11974 IChevaliej[l977l : iMcCrav fc Snowlll979l : ISpitzeiiliggof T 
Canonically, only a tenth of a typical supernova's kinetic 
explosion ener gy is available for feedback, since the rest is 
radiated away (|Larsonlll974 . but the debate about the un- 
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of the details in this paper do rely on models which treat 
ULXs as stellar-mass BHs. We mostly consider populations 
of XBs rather than individual sources, and in young stel- 
lar populations then XBs with more massive donor stars 
(including, we expect, ULXs) should dominate the luminos- 
ity of the XB population; most of this work discusses the 
potential effects of feedback from these sources. In old stel- 
lar populations then LMXBs are still present, and we also 
briefly consider to what extent LMXBs might keep gas hot 
in non-star forming galaxies. 

The importance of binary stars on the appear- 
ance of stellar populations, and therefore on the 
appearance of gala xies, has been demonstrated be- 
fore 
19971: 



( see, e.g., Vanbeveren. van Bever fc De Dondeij 
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derlying physics and characteristics of SNF is ongoing. For 
example, ICeverino fc KlvpinI l|2009l 'l argued that SN feed- 
back may be dominated by the supernovae from runaway 
stars which are able to explode in cold gas, since SNe ex- 
ploding in hot gas lose more of their energy to radiation. 
The fact that around 40 per cent of luminous core-collapse 
supernovae se em to occur in runaway s tars, as deduced from 
observ ations llJames fc And erson 2006) and population syn- 
thesis (jEldridge. Langer fc Tout 2011 ) supports the poten- 
tial significance of this effect. 

Assu ming that SNF is dominant, the classic work of 
iDekel fc Silk (198^) found that 'The critical condition for 
global gas loss as a result of the first burst of star for- 
mation is that the virial velocity be below a critical value 

on the order of 100 km s~^'. However, some of the lowest- "t~. ":~i — . - .'. . 

mass dwarfs (with virial velocities ^10 km s~^) appear to iBelkus et al.' '2003'; Zhang et al. l2005l: iDionne fc Robert! 

somehow r etain gas after th e first episode of star formation 
(see, e.g., Ilrwin et al.l 120071 ): hence in some cases feedback 
seems remarkably weak. Yet feedback may often need to be 
strong in order to match the observed dwarf galaxy popula- 
tion with the predictions of CDM (see, e.g.. IKoposov et al.l 
[2009) . Continuing inadequacies in our understanding of stel- 
lar feedback i s one o f the possible explanations suggested by 
ISawala et al.l (|201in for the fact that models of galaxy for- 
mation and evolution in « 10^'' M© halos produce results 
apparently inconsistent with observations. For such reasons 
we must continue to examine all possible feedback mecha- 
nisms. 

In this paper, we argue that another source of stel- 
lar feedback - X-ray binaries (XBs) - should be included 
in galaxy modelling (see also ICantalupol l20ld ). X-ray bi- 
nary feedback (XBF) should be especially important at low 
metallicity and hence in the early e pochs of galaxy form ation 
(see also the independent work of iMirabel et al.ll201ll ). We 
emphasise how the properties of XBF seem qualitatively dif- 
ferent to those of SN feedback, and that these differences are 
interesting and important. We also argue that the stochastic 
variation in XBF should be expected to be more significant 
than for SNF. Partly for that reason, XBF seems promising 
in helping to explain the complex star-formation histories 
and diversity of very low-mass dwarf galaxies. We also note 
that strong XBF can begin before strong SNF, which may 
allow XBF to modify the effects of SNF. 

The class of XBs contains objects with a wide range 
of ages and luminosities. Accreting stellar-mass black 
holes (BHs) and neutron stars were once broadly di- 
vided into two groups; low-mass X-ray binaries (LMXBs) 
where mass is transferred from the low-mass (< 2 M© 
donor star) by Roche-lobe overflow and high-mass X- 
ray binaries (HMXBs) where the high-mass (>10 M0) 
donor star spontaneously loses mass in a wind or ex- 
cretion disc. The majority of luminous systems known 
as ultra-luminous X-ray sources (ULXs) are probably 
XBs, but in neither of those old categories; they are 
expected to contain intermediate-mass or massive donor 
stars which are under going short-lived phases of Roche- 
lobe overflow (see, e.g., lGilfanov, Grimm fc Su nvaevll2004al : 
iRappaport. Podsiadlowski fc Pfahl 200 5 , hereafter referred 
to as RPP; iMadhusudhan et all I2OO6I : iBlecha et all I2OO6I : 
IZezas et al.ll2007l ). The nature of ULXs does not affect our 
main conclusions, which are consistent with empirical stud- 
ies of XB populations, as we describe later. However, some 
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lEldridge fc Stanwavil2009i V The potential effect on galactic 
structure of a subset of compact binari es other than XBs 
has a lso recently been considered by iMoore fc Bildstenl 
l|201ll ): binary stars matter on galactic scales. 

Section [2] outlines some of the reasons why we con- 
sider that energy input from XBF could be important in 
shaping galaxies. This includes a brief summary of the ob- 
served effects of some XBs on the ISM, the likely metal- 
licity dependence of the formation rate of young XBs, and 
how XBF can operate before SNF. There we also outline 
why we think that XBF may be particularly relevant to 
some of the observed characteristics of dwarf galaxies. Sec- 
tion [3] expresses in more detail the energy available from 
mean populations of XBs, and provides expressions for use 
in semi-analytic galaxy models. In Section|3]we examine po- 
tential stochastic variations in XBF, and suggest that they 
provide a natural reason for diversity in the evolution of 
initially identical dwarf galaxies. In particular. Section 14.11 
flnds that the most important factor in determining the ef- 
fect of XBF on individual low- mass halos may not be metal- 
licity but probability. Furthermore, in Section [S] we argue 
that natural statistical variation in XBF may explain why 
the effect of XBF would be most commonly significant in 
populations of galaxies with stellar mass approximately be- 
tween 10^ and 10* M©, potentially resolving the problems 
wit h current gal a xy ev olution models which were identified 
by Sawala et ajj (|201l|). Our discussion speculates on how 
the interaction between XBF and SNF may lead to a rich 
variety of behaviour. 



2 MOTIVATING X-RAY BINARY FEEDBACK 

2.1 Efflciency of energy coupling & demonstrable 
outflows 

A large fraction of the energy output of at least some XBs 
seems to be deposited into the ISM. By using ionisation neb - 
ulae around ULXs as calorimeters, IPakull fc Mirionil (|200J) 
found that the energy deposited in the nebulae was compara- 
ble to the observed X-ray luminosity ( assuming isotropic X - 
ray emission; see also ULX IC 342 X- 1. [Roberts et al.ll20o3 '). 
The potential photoionisation effect of XBs on galactic- 
scale stella r feedb ack has been indep e ndent ly considered by 
ICantalupol (|2010l ') and lMirabel et al.l l|201lh . 
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Photoionisation is not the only potential route through 
which XBs might affect their environment. Other XBs also 
seem efficient i n depositing energy in the ISM to produce ki- 
netic outflows: iGallo et al.l ([2005) found that 'the total power 
dissipated by the jets of Cygnus X-1 in the form of kinetic en- 
ergy can be as high as the bolometric X-ray luminosity of the 
system' . The nebula around SS433 (W50) ap pears to be jet- 
inffated, and SS433 drives outflows (see, e.g., Blundell et al.l 



I2OOII : iGoodafl. Alouani-Bibi fc BlundeilbOllh. The size of 



the X B-inflated nebula observed bv lPakuU. Soria fc MotchI 
l|2010l) is ~ 300 pc. This kinetic input from XBs wa s con- 
sidered by iFender. Maccarone fc van KesterenI (J2005l l as a 
potential driver for cosmic-ray emission on galaxtic scales. 

Indeed, some XBs put far more energy into 
local outflows thaii into their X-ray luminosity: 
IPakuU. Soria fc MotchI (|2010t ) argue that the kinetic 
output of one microquasar is ~ 10** times its X-ray power. 
This XB is not a ULX, but its kinetic energy ou tput is a 
few times 10*° erg s~^ (see also ISoria et al.ll201Cl ). If such 
efficient deposition of energy into the ISM is common, then 
we vastly underestimate the magnitude of XBF in what 
follows. In that case then the majority of the feedback 
would be in kinetic outflows rather than photoionisation. 

2.2 Metallicity dependence 

The abundance of bright, young BH XBs (in particu- 
lar ULXs) is exp e cted to increase at lower metallicity 
(Z). iLinden et al.l (|2010l ) and Justham & Podsiadlowski 
(in prep.) have studied the metallicity-dependence of lu- 
minous BH XBs, roughly agreeing that the formation rate 
at Z ~ 0.01 Zq could be ~10 times higher than at Zq. 
Explaining the detailed reasons for this metallicity depen- 
dence is left to those papers, but it is not enough to ar- 
gue that BH formation from single stars should be more 
common at low metallicity; the changing formation rate of 
interacting binaries containing BHs as a function of metal- 
licity turns out to be dominated by changes in the available 
range of initial binary separations which can go on to pro- 
duce such systems. Specifically, the main effect relates to 
the allowed timing of the interaction when the envelope of 
the BH-forming star is removed (with respect to the nuclear 
evolu t ion o f that star). For non-conservative mass transfer, 
iDravl Hooi) found a BH HMXB formation rate at Z = 0.001 
that is ~10 times higher than ai Z — 0.02. Earlier work by 
iBelczvnski. Sadowski fc Rasid (|200J) had also found that 
the number of young binaries consisting of a BH with a 
main-sequence star (BH+MS) should change by a large fac- 
tor as a function of metallicity. Belczynski et al. looked at 
model populations with metallicities oi Z = 10~*, 10"'' and 
0.02. At an age of 11 Myr, the ratio of BH-I-MS binaries 
between these models was « 9 : 5.7 : 1, and at 104 Myr, ~ 
6.7 : 4.5 : 1. The production rate of BH-I-MS systems is not 
the same as that of ULXs, but it should be indicative of the 
production rate of bright BH XBs. 

All the XB formation calculations discussed in this sub- 
section have assumed that the low-metallicity IMF is the 
same as the present-day IMF, so a top-heavy IMF at low 
metallicity would enhance this metallicity dependence. (We 
note in passing that binaries themselves should affect ob- 
serva tional determinations of the IMF: see, e.g., lEldrided 
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Figure 1. The black curves show estimates of the time available 
between BH formation and the start of the most energetic SNF, 
as a function of the initial stellar mass which produces the BHs, 
and assuming that SNF arises from stars below 2OM0 (see text). 
During this delay XBF could operate before SNF. The grey curves 
show the approximate lifetimes of massive stars up to 4OAf0 , from 
which the black curves are calculated. The horizontal dotted grey 
line represents the (shorter) delay between star formation and 
ULX onset found by RPP, which indicates that we may have 
underestimated the time gap available for XBF before SNF. 



Furthermore, the typical masses of stellar BHs m ay well 
be h igher at lower metallicity (see, e..g., [Eldridgc fc Touti 
I2OOJ). Any increase in the typical mass of accreting BHs 
at lower metallicity should increase the pote ntial luminos- 
ity o f the XB population further (see, e.g., iLinden et al.l 
[2OI0I). We note that a novel formation mechanism has re- 
cently been proposed for a subs et of XBs containin g un- 
usually massive stellar-mass BHs (|de Mink et al.ll2009l ), but 
any metallicity dependence of the rate for this alternative 
format ion channel is so far u nclear. From observational evi- 
dence, ICrowther et al.l (|2010l ) point out that the BH masses 
in known wind-accreting BH HMXBs do seem to be larger 
at lower metallicity. 

We are not aware of an observed sample of ULXs with 
population II metallicity, so the predicted range of ULX 
metallicity dependence is currently hard to test. However, 
IZampieri fc Roberts! ((20091) concluded that even in the avail- 
able local samples of UL Xs there are hints of a metallic- 
ity dependence (see also iMapelli. Col pi fc Zampieri 2009|: 
Mapelli et al. 2010). In addition. IShtv kovskiv & Gilfano^ 
(200f |) compared the small ma g ellanic cloud (SMC) to 
the .Grimm. Gilfanov fc SunvaevI (|2003r ) relation between 
HMXB populations and the SFR, as derived for more mas- 
sive, higher-metallicity galaxies. Using standard far-IR, Ha 
and UV indicators of the SFR, the HMXB population in 
the SMC appears be overabundant by a factor of ~10 (see 
also Dray 2006). If that discrepancy is due to a metallicity 
dependence, then it is significantly larger than the one we 
apply in this paper. 
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2.3 The first supernova may not be the first 
energetic supernova: quiet black hole 
formation 

In a given starburst, BHXBs may be able to provide en- 
ergy input before the first energetic supernova occurs. From 
slightly below solar metallicity and lower, the majority of 
single stars more massive than « 40 Mq are expected to 
produce black holes by direct collapse, and the associated 
SNe are expected to produce no significant explosion, per- 
haps contributing no energy at all to f eedback (see, e.g. , 
Frvej Il999l: iMirabel fc RodriguesI l2003l : iHeger et al.l booj : 



Eldridge fc ToutI |2004| '1. SNe which produce BHs by fall- 



back may also be very faint and contribute little to super- 
nova feedback (a pos sible example of su ch a supernova is 
1997D, for which see iTuratto et aLlliggsT). A t a metalhcity 
of ^ 0.001, figure 5 of lEldridge fc ToutI (|200i ) indicates that 
direct collapse occurs for a ZAMS mass of ^ 35 Mq , and BH 
formation by fallback for masses ^ 20 M0 . This also means 
that the normalisation of the energetic SN rate to the SFR 
and that of the total SN rate to the SFR are not identical. 
So most stars ^ 20 Mq may hardly contribute to feedback. 
Of course, the previous description of events which con- 
tribute to SNF is simplified. Fir stly, binary evolution makes 
the si t uation more complex (e.g. Podsiadlowski. Joss fc Hsul 
1992"; 'Do Bonder & Vanbevercn' '2003'; 'Eldridge^ l2005l : 
Eldridge, Izzard & Tout 2008). In addition, extremely mas- 
sive stars are expected to pro duce pair-instability supernova , 
especially at low metal l icity (IBarkat. Rakavv fc Sackl 119671 : 
iRakavv fc Shavivlll967l : iHeger et al.ll2003l '). and at least one 
convinc ing candidate for su ch an explosion has been ob- 
served (jGal-Yam et al.ll2009l ) . Recent examples of unexpect- 
edly luminous SNe also indicate that fur ther exceptions to 
the simple pictur e above probably ex ist fSmith et al."2007; 
Ofek et all l2007l: [Quimbv et al 
Gezari et al.' '20091: Miller et al 
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Rest et al. 2011 : Quimbv et al.ll201ll ) . These have often been 



linked with luminous blue variable progenitors, i.e. stars sig- 
nificantly more massive than 20 Mq, although that is not 
consistent with prior expectations. We note that the lead- 
ing explanation for the high observed luminosity of most 
of these events involves a SN which is more efficient at ra- 
diating awa y kinetic energy, not an intrinsically more en- 
ergetic SN jSmith fc McCravl l2007l : lAgnoletto et^ l2009l : 
ISmith et al.ll2010|). Those excep t ionall y luminous events are 
rare, although iGal-Yam et al.l (|2007l ) suggest that many 
more events - all Type Iln SNe, i.e. ~ 2 per cent of all bright 
core-collapse SNe - might be from such massive stars. It is 
possible that a non-negligible fraction of SNF begins sooner 
after a starburst than suggested by the arguments in the 
previous paragraph. 

Despite those caveats, we consider it is very plausible 
that, in the majority of starbursts, the first BHs could be 
formed well before many 'loud' supernovae from that star- 
burst explode and drive significant SNF, especially at low 
metallicity and hence in the early universe. In principle, the 
first BHXB could switch on almost as soon as the first of 
those BH is born. The limit is the evolutionary timescale of 
the donor star: mass transfer should be stable from a radia- 
tive donor up to a few times more massive than the mass of 
the BH. Figure [1] indicates the approximate amount of time 
which might be available to XBF before strong SNF begins. 



Table 1. Approximate supernova feedback efficiency necessary 
for starburst gas ejection. 

VvMai \u\/j:esn 

(km S-l) A/tot/Mstollar=5 Aftot/Mstcllar = 50 



100 
10 

1 



0.1 
10-3 
10^5 



1 

0.01 
10-* 



Stellar lifetimes are e stima ted using equations 4 and 57 of 
iHurlev. Pols fc Tout! (|2000l '). This suggests that if SNe from 
progenitors more massive than 20 Mq contribute little to 
feedback then XBF should have several Myr of free time to 
act before SNF begins. If the time of ULX onset found by 
RPP is correct, then there is almost 6 Myr left before 20 Mq 
progenitors explode and «1 Myr before 40 Mq progenitors 
do. During this time XBF could modify the environment on 
which SNF acts. 

Note that, in adopting the work of RPP, we have implic- 
itly adopted their assumption that BHs are born with a mass 
of 10 Mq, along with their limited range of donor masses; 
we expect that the time before the first ULX switches on is 
likely to be even shorter than the minimum time given by 
RPP. Moreover, if stellar BH masses increase with decreas- 
ing metallicity then the maximum stable donor star masses 
also increase, which could further shorten the time before 
XBF can begin. 



2.4 Gas retention in dwarf galaxies 

Feedback in some dwarf galaxies seems to be able to 
truncate SF without ejecting all the gas, allowing ex- 
tended or repeated SF episodes even in some astonish- 
ingly tenuous dwarf galaxies (e.g. Leo T, for which see 
Irwin et al.l l2007l: iRvan-Weber et al.l boosi : Ide Jong et all 
2008J: iGrcevich fc Putmanll2009l ). Man y dwarf galaxies do 
lose their gas (for a review see, e.g., iTolstov. Hill fc Tosil 
12009!) . but the evidence that some do not is much more sur- 
prising, since for galaxies with virial velocities <10 km s^^ 
even very low SNF efficiencies could suffice to unbind the 
gas. We illustrate this in Table [1] in which we give order-of 
magnitude estimates for the efficiency of SNF which would 
be necessary to eject all the gas from low-mass starburst 
galaxies. We show estimates for two different ratios of total 
mass (Mtot) to stellar mass (Afstoiiar), one of which repre- 
sents an extremely dark-matter dominated case. We simply 
relate the binding energy (|L'^|) to the virial velocity (Vviriai), 
taking \U\ = Mtot Kdriai- Our estimates also assume that one 
SN contributes to SNF for every 100 Mq of stars formed 
in the starburst, and that perfect SNF efficiency provides 



10^ erg/SN. The important point is that SNF could be less 
than 1 per cent efficient and would still eject the gas from 
realistic dwarfs with virial velocity < 10 km s-^. 

Energy injected on the long timescale of XBF input 
(~10* yr) seems more likely to be able to heat the gas reser- 
voir non-destructively than energy injected rapidly com- 
pared to the dynamical timescale of the galaxy (as in SNF) . 
Since we also expect that XBF can begin before strong 
SNF, we furthermore suggest that XBF might create 'chim- 
neys' to help channel subsequent SNF away. Previous stud- 
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ies of how SN explosions eject or enrich gas reservoirs from 
dwarf galaxies have sometimes found it possible that the hot, 
metal-enriched gas can escape through such chimneys whilst 
leaving more of the col d phase intact than simpli stically ex- 
pected (see especially [Mac Low fc Ferrara 19991: also, e.g. . 



be Young fc Gallaghed 119901 : iDe Young fc HeckmanI il994l : 
Fragile et aLlboOSl ). Such previous work has found that SN 
remnants preferentially expand through regions containing 
hot, low-density gas than through the colder, higher-density 
phase outside this 'chimney', and the SN ejecta can thereby 
be channelled out of the galaxy. Instead of a situation where 
early SNF creates chimneys for later SNe (as in previous 
work), XBF may be able to form such hot, low-density re- 
gions before significant SNF begins. Even if XBF simply pre- 
heats the ISM before SNF, without helping to form chim- 
neys, the n it might lead to a red uction in SNF efficiency 
(see, e.g.. lCeverino fc Klvpinll2009l . and references therein). 
Arguments in favour of a relativ ely gentle f e edbac k pro- 
cess in dwarf galaxies were made bv lLee et al.l (|200q ). who 
suggested that 'a less energetic form of metal-enhanced mass 
loss than blowouts could explain the small scatter in the 
L — Z and M« — Z relationships'; XBF may help fulfil that 
need. Furthermore, the stochasticity of XBF, as presented 
in Section 3] may help explain the diversity in the known 
Milky Way satelite galaxies, e.g. how Leo T can differ so 
much from the others in present-day neutral gas content. 

2.5 Timescales Sz episodic star formation 

A mechanism which stops star formation (SF) without 
ejecting the gas reservoir may allow episodic SF. To ex- 
plain separated bursts of SF seen in some dwarf galaxies, 
iQuillen fc Bland-HawthornI l|200a ) suggest that the feedback 
in those galaxies should be delayed over a timescale of 50- 
100 Myr. The delay timescale is not the only important fac- 
tor in their model, but they find that such a delay is nec- 
essary to produce episodic star formation. As there seemed 
to be no stellar source of feedback with a suitable timescale, 
they concluded that the explanation for this timescale must 
be non-stellar. Based on the models of RPP, energy input 
from XBF can increase for ~10* years after the onset of SF, 
apparently providing a stellar source of feedback naturally 
consistent with that model for episodic SF (see also Section 
[Sjl. Continuing energy input from longer-lived XBs might 
further delay the subsequent star formation episodes. 

Stellar feedback could perhaps be delayed in another 
way such that it could satisfy the conditions found by 
iQuillen fc Bland-HawthornI l|2008r ): they themselves specu- 
late on how the effective timescale of SNF might be ex- 
tended. Some of the feedback from SN la could happen 
on a sufficiently long timescale and some core-collapse SNe 
could be delayed b y the effect of binary evolut ion (see, e.g., 
De Ponder fc Vanb cv crcn 2003; Eldridgc, Lane:er fc Toud 
2011^ . The model of IQuillen fc Bland-HawthornI l|2008l )"ls 



not the only attempt to explain ho w episodic SF can oc- 
cur. F or example, the simulations of iPontzen fc Governatd 
l|2011J) produce bursty star formation in a dwarf galaxy with- 
out modifyi ng their input from stellar feedback. Nonetheless, 
based on the lQuillen fc Bland-HawthornI (|2008h model, XBF 
may help to explain this effect; we also repeat that it seems 
hard for simple SNF to allow gas retention in galaxies such 
as Leo T. 



The extended duration of energy input from XBF also 
means that XBs born with low metallicity in a galaxy with 
a large specific SF rate could still be active when the star- 
forming metallicity of the galaxy has already been enriched; 
this could be important if the XBF input is strongly metal- 
licity dependent. 



3 MEAN X-RAY BINARY POPULATIONS 

Here we present the expected energy output from mean X- 
ray binary populations, as measured by their X-ray lumi- 
nosities. We have already argued that the efficiency with 
which these outputs might couple to the ISM could be very 
high (see Section [2. ip . The potential importance for galaxy 
evolution of the mean energy inputs described here are illus- 
trated by scaling relations in the appendix. However, much 
of the later sections discuss how these simple treatments 
alone may not capture the essence of XBF, especially for 
low mass or low star formation rate (SFR) populations. 



3.1 X-ray binaries in young populations 

iGrimm. Gilfanov fc SunvaevI (|2003l ) observationally found 
that the most luminous sources dominate the X-ray out- 
put of a star-forming galaxy, so here we can concentrate 
on those most luminous sources when describing the power 
output from young XB populations. We use the theoreti- 
cal models of RPP, which assume that those dominant X- 
ray sources are ULXs formed from stellar-mass BHs, but 
we also demonstrate that those are consistent with observa- 
tions. In Section|3|we consider the situation at low absolute 
SFR, when it is unlikely that the most luminous source in a 
galaxy is a ULX. 

Scaled to a steady-state SFR of 1 core-collapse SN 
per century, RPP's models produce a total accretion lu- 
minosity from ULXs, (Lulxs) of 2.3 x lO'^^erg s^^ if the 
Eddington luminosity can be exceeded by a factor of 10, 
and 1.5 x lO^^erg s~^ if the Eddington luminosity is strictly 
obeyed. Their calculation was performed at solar metallicity 
(Zq), and assumed BHs of mass 10 M0. Whilst these cal- 
culations are subject to several uncertainties, RPP seemed 
to reproduce well the observed ULX distribution in the 
Cartwheel galaxy. Based on their results at Zq we assume: 



LxBF « 2 X 10^" erg s'^ 



SFR 

M© yr- 



(1) 



for the luminosity of the entire XB population in a mean 
star-forming galaxy. This seems consistent with the relation- 
ship between SFR and total X- ray luminosity in the 2-10 keV 
band (Lx,2-iokov) found by ICrimm. Gilfanov fc SunvaevI 
(2003). In the linear regime (i.e. for mean populations), their 
empirical relation is: 



ix,2-i0kcv = 6.7 X lO^'' erg s ' 



SFR 

Meyr- 



(2) 



which would be equivalent to equation [Tj if « 2/3 of the 
bolometric luminosity is not observed in the 2-10 keV band. 
It also seems likely that the SFR normalisations are not 
identical when comparing the empirical to the theoretical 
luminosity. 
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Figure 2. The three black curves show the energy output from the young XB population of three different SFHs, i.e. an instantaneous 
burst of SF, a constant SFR, and an exponentially decaying SFR. Each simple population contains a total SF mass of IO^Mq, at constant 
metallicity (Zq). For comparison, we also show grey lines showing the SN energy input (assuming one SN per 100 Mq of SF, and 10^^ erg 
of energy per SN), but for 'ongoing', i.e. pre-existing, constant SF rather than for the new SF which produced these XB populations; for 
a discussion of the likely delay to SNF, see section [2]3] and Figurc[l] The standard SNF line assumes the canonical 10 per cent efficiency. 



We take Equation [T] as potentially giving the energy in- 
put to XBF at Zq. This is somewhat less than the approx- 
imate energy output of one SN (~ 10^^ erg) per century 
(fs 3 X 10^ s), i.e. « 3 X 10*^ erg s~^, which is roughly ap- 
propriate for a SFR of 1 Mq yr~^ (depending on the IMF). 
However, XBF would become competitive if the efficiency 
of XBF is greater than of SNF (as previously suggested). 
Furthermore, since we expect that the importance of ULXs 
increases as the metallicity decreases (see Section [2.2^ then 
ULXs could easily dominate SNe as an energy source for 
feedback as the metallicity decreases. 



3.2 X-ray binaries in old populations 

LMXBs may also provide u seful energy inpu t in old popu- 
lations to keep gas hot. The iGilfanovl ( 2004a ) relation: 



Llmxb ~ 10 erg s 



f M^ 



VlQS Mq 



(3) 



implies significant energy input in a massive elliptical (see 
the scaling relations in the appendix for estimates of how 
much gas might be kept hot by XBF in old populations). 

The metallicity dependence of this energy input is 
less clear than for ULXs, though there is no expecta- 
tion that it should be as strong as for young BHXBs. 
For globular clusters (GCs), XB luminosity is higher 
in metal-rich GCs than metal-poor one s, but the ori- 
gin o f this effect is unsolved (se e, e.g., iBellazzini et al 



19951 : 

20041: llvanoval l2006l : 

Sivakoff et al.ll2007l ). 



Kundu. Maccarone fc Zepj I2OO2I: Ijordan et al 



Kundu. Maccarone fc Zepj l2007l: 



3.3 Mean X-ray binary feedback prescriptions for 
semi-analytic models 

We now suggest a simple treatment of XBF intended for use 
in semi-analytic models, based on our earlier descriptions. 
Here we assume that XBF efficiency (measured with respect 
to XB luminosity; see Section l^TII is 100 per cent, but clearly 
another parameter could be added. 

The time- and metallicity-dependence of the energy in- 
put from young XBs makes their inclusion slightly compli- 
cated. We have produced a fit to the energy output from a 
ULX population produced by an instantaneous burst of star 
formation (using model B in figure 8 of RPP). This fit is in 
two parts. We use the fitting function I/xbf,i for the XBF 
energy input at early times, when there is a high-luminosity 
plateau in the power output. At later times, we use I/xbf,2 
to represent the decline of the XB power output as the pop- 
ulation ages. For even older stellar populations, longer-lived 
XBs should set a floor on the XBF energy input. For mean 
populations. Equation |3] could be used to give the minimum 
energy input from XBs at any age. 

For any star formation history (SFH) other than an in- 
stantaneous burst then at some time tnow the instantaneous 
power input from XBF will be a combination of contribu- 
tions from populations with different ages. For young pop- 
ulations (< 1 Gyr) the XB power output should be the sum 
of both LxBP 1 and I/xbf 2, where: 



L^ 



erg s" 



2 X lO'* 



SFR{t') 



1 M 



yr" 



X f(Z{t'))dt' (4) 



and SFR(t') and Z(t') were the star formation rate and 
metallicity of the newly-forming stars at some previous time 
t'. We require /(Zq) = 1.0, and suggest setting /(Zq/200) = 
10. Note, however, that this function is uncertain; further 
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observational and theoretical work is necessary to determine 
it more accurately. 

We estimate that ti — 10 ' and i2 = 10 ' years (again 
using model B in figure 8 of RPP). As explained in Section 
12.31 this value for ii could probably be reduced. At later 
times we approximate the power input from this young com- 
ponent as steadily declining: 



JXBF,2 _ , „57.85 



erg s" 



(tn 



t'Y 



SFR{t') 



1 M, 



-dt' 



yr- 



(5) 

(we have suppressed the metallicity dependence here for 
brevity, but presume that the same function as used in Equa- 
tion [4] could be applied here). 

Figure [2] shows the outcome of the above equations at 
solar metallicity. We re-emphasize that, at low metallicity, 
the relative importance of XBF in young populations should 
increase. As a result then, if the self-enrichment time of a 
particular galaxy is shorter than ~ 10* years, then the to- 
tal feedback from XBF might easily still be increasing even 
when the mean metallicity of the stars in the galaxy is no 
longer low. 



4 THE IMPORTANCE OF STOCHASTICITY 

Galaxies with low absolute mass and SFR are likely 
to deviate from the mean relations given in the pre- 
vious section. The effect that random sampling of the 
initial mass function has on observed stellar popu- 
lations has been in ves tigate d by many authors (see, 
e.g., ICervino et aD I2OO2I: ICervino fc Valls-Gabaudl 
20031: IVillaverde. Cervino fc Luridianal I2OIC 



da Silva. Fumaealli fc Krumhold I2OI2I : lEldridg"3 mu 
When considering the formation of XBs, the parameter 
space of binary initial conditions is multi-dimensional, and 
the volume of that space which produces bright XBs is 
much smaller than the volume which leads to a massive 
star and a supernova. So we expect stochastic variation to 
be greater for XBF than for SNF, especially because the 
power output from XBs should be dominated by a small 
fraction of the most luminous sources. Here we estimate 
the potential importance of this effect, in particular by 
comparing the distribution of Lx, the XBF luminosities for 
each population within an ensemble, to <Lx>, the global 
mean value. 

At the peak luminosity of the population synthesised 
by RPP, the number of luminous BHXB systems in a star- 
burst which produces 10^ SNe is «30, whilst at an age 
of 10^ yr that population only contains ~10 sources with 
L > 10'^^ erg s~^ (there is a small model-to- model variation 
between the numbers of expected systems, as shown in their 
figures 6 & 7, but not one which is large enough to affect this 
argument). Because the majority of the expected XBF lumi- 
nosity at those early times is produced by far fewer sources 
than lead to the SNF energy input, statistical variation be- 
tween populations at these times should be more significant 
for XB F than for SNF. 

iGilfanov. Grimm fc SunvaevI (J2004bf l showed more for- 
mally that the variation between young XB populations in- 
creases significantly as the SFR decreases. They analytically 
combined the probability distribution function for the num- 
ber of sources present with the luminosity function from 
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Figure 3. The top panel shows how the statistical variation in 
XBF populations may be realised for diflerent absolute SFRs (at 
Zq); see text for details. The lower panel uses the cumulative 
distribution of the same population models to show the increased 
importance of the /lig/i-luminosity tail of the distribution for low 
SFRs. As the mode of the distribution falls for low SFR, the 
effective range of XBF input in low-SFR halos extends over two 
orders of magnitude, even ignoring metallicity effects. 



which those sources are drawn; they thus found the proba- 
bility distribution for the total luminosity (below we demon- 
strate this numerically). They found that the regime in 
which small-number statistics affects the expected total lu- 
minosity extends much further than might naively have been 
expected; e.g. for populations of young X-ray binaries it can 
extend to SFRs of a few solar masses per year. This may 
well have implications for the evolution of populations of 
low-mass galaxies, in particular since the probability distri- 
bution functions derived by Gilfanov et al. become strongly 
asymmetric at low SFR. 

4.1 Low star- formation rate effects for young 
populations 

For the reasons given above, it seems possible that stochastic 
differences in XBF magnitude between galaxies with simi- 
lar initial properties - i.e. random halo-to-halo variations 
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- might be important when trying to understand the ef- 
fects of XBF. Here we estimate the potential significance 
of XBF stochasticity in young populations by producing a 
probability distribution function for the XBF magnitude. 
To do this we use a Monte-Carlo rea lisation of the work of 
iGilfanov. Grimm fc SunvaevI (|2004a ). i.e. we populate each 
simulated starburst with a collection of XBs which are ran- 
domly selected from an X-ray luminosity function (XLF) in 
such a way as to preserve the appropriate mean for the total 
luminosity in the large-SFR limit. We adopt their power-law 
XLF with an index of —1.6, and assume that no ULX can be 
more luminous than 10^^ erg s~^. We further allow the XLF 
of the individual sources to extend over eight orders of mag- 
nitude for these young populations; this choice affects the 
details of the results, but not the general principle. Specifi- 
cally, reducing this range over which the power-law operates 
increases the fraction of galaxies with no XBF whatsoever, 
rather than with little XBF; that makes for uglier plots, but 
does not change the qualitative conclusions. Figure |3] shows 
the outcome of those calculations, demonstrating that the 
bulk of halos with low absolute SFR could put out X-ray 
luminosities an order of magnitude below the mean value, 
<Lx>, whilst several percent of initially similar halos could 
experience XBF an order of magnitude stronger than the 
mean value. This two-order-of-magnitude spread is larger 
than the systematic increase in mean XBF magnitude, by a 
factor of roughly 10, which we expect to occur at low metal- 
licity. 

For low-mass, low-absolute-SFR galaxies at low metal- 
licity, the effects of metallicity and stochasticity should com- 
bine, meaning that for a small percentage of halos the spe- 
cific XBF (relative to the SFR) could be two orders of mag- 
nitude higher than the mean value given by Equation [T] 
Note that as the expected number of bright sources increases 
within the galaxy as the metallicity drops then the spread at 
a given SFR also decreases; ignoring any contribution from 
higher mean BH masses then the shape of the Lx distri- 
butions is constant for total <Lx>. It is possible that the 
shape of the XB luminosity function is also a function of 
metallicity; if so this may be as important as the metallicity 
dependence of the XBF normalisation. 

In high-absolute-SFR populations, the likelihood of a 
large deviation from the mean is small, and the mean re- 
lations from section [3] can be applied. The importance of 
stochasticity to XBF should also be somewhat reduced for 
situations in which the significant input from XBF can be 
averaged over a relatively long time, i.e. over several genera- 
tions of luminous XBs. However, when trying to understand 
the diversity in the baryonic content of low-mass halos, we 
suggest that probability distributions for XBF, such as that 
described above, should be adopted. We emphasize that this 
should be true even for high specific SFRs. 



4.2 Low mass effects for old populations 

We do not repeat the above calculations for LMXBs, as 
the remainder of this paper will be mostly concerned with 
XBF in young populations. However, we will briefly con- 
sider when XBF in old populations might be affected by 
significant stochastisicity. Assuming a quantised LMXB lu- 
minosity of 10^^ erg s~^. Equation [3] gives a rough expec- 
tation value of one LMXB for a mass of lO^Mo. The more 
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Figure 4. We estimate the relative importance of XBF for modal 
galaxies in the XBF probability distribution by taking the ratio 
of the modal XBF luminosity (I/x) to the gravitational binding 
energy ("l?Galaxy) for each galaxy mass. We expect that the low- 
Z curve (broken black line) should be relevant for the galaxy 
mass range at which this ratio peaks. As this figure only indicates 
the significance of XBF for the mode of the galaxy population, 
we repeat that the shape of the XBF probability distribution 
itself is relatively narrow and roughly symmetric for galaxies more 
massive than those at the peak in this figure, but increasingly 
broad and asymmetric for less massive galaxies (see Section 14.11 
and Figure 1311. 



sophisticated analysis of lGilfanovl (|2004al ) finds an expected 
number of LMXBs with luminosity > 10^ ^ erg s~^ oi ^ 140 
per 1O^^M0. Based on iGilfanovl (1200481 ) ■ we suggest that 
using mean values of LMXB XBF might be misleading for 
galaxy masses nruch less than IO^^Mq, so a PDF for 'old' 
XBF may also be needed to faithfully reproduce systems 
where LMXBs are significant. 

As canonical LMXB lifetimes are ~1 Gyr, then any 
dwarf galaxies low enough in mass to posess (typically) one 
significant LMXB could be expected to experience heating 
episodes lasting from a few hundred Myr to several Gyr. 
Cooling and a new epoch of star formation might be de- 
layed until the LMXB heating stops. Scaling relations in 
the appendix suggest that it may be reasonable for a bary- 
onic gas fraction of several tens of percent to be kept hot 
in some low-mass galaxies by XBF. The inherent random- 
ness of XBF may again be important in understanding the 
present-day diversity of dwarf galaxies, and should be con- 
sidered in future models. 



5 FOR WHICH GALAXY MASSES SHOULD 
X-RAY BINARY FEEDBACK BE MOST 
COMMONLY SIGNIFICANT? 

As explained in the previous section, then below a certain 
absolute SFR the most probable XBF input for a galaxy 
with a given SFR should be below the expectation value 
for XBF energy input at that SFR. That is, the mode of 
the XBF probability distribution is below the mean. In 
that low-SFR regime, the location of the mode scales as 
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(SFR)^/^ for an XLF slope of -1 6 (see Sectfon O and 
iGilfanov. Grimm fc Sunvaevll2004a ). Above this regime, i.e. 
for high enough SFRs that the mode of the I/x/<ix> dis- 
tribution is very close to <Lx>, then clearly the most likely 
XBF input scales linearly with the SFR. This allows us to 
make a simple estimate of which part of the galaxy popula- 
tion might be the most likely to be modified by XBF. 

We will assume that th e gravitational binding energy 
typic ally scales as M^"'^ (e.g. lSaitolll979l : lYoshii fc Arimotol 
119871 ). and that we can take a roughly constant star forma- 
tion lifetime AisF. In this case then the galaxy mass where 
the ratio of XBF energy input over gravitational binding 
energy is highest - for a typical halo, rather than for the 
tail of the distribution - will arise at the transition between 
these low-SFR and linear regimes. We can use these assump- 
tions to estimate the ratio between the XBF luminosity and 
galaxy binding energy for 'typical' galaxies of different star- 
burst masses. Figure 1X21 shows these estimates, assuming a 
single dominant starburst with a duration of 10* years. 

Using our previous assumptions, and writing the effec- 
tive upper end of the XLF as Lmax then the stellar mass M* 
corresponding to the point of most common, most significant 
XBF is estimated to be: 



M, 



^^xWMqJ VlxlO«ergs-i 
at solar metallicity, or: 



5 X IO6M0 



1 X lO^Oerg s- 



A tsF 

4 X 10'' yr 



4 X WyrJ 



(6) 



(7) 



if we allow a factor of ten enhancement to the number of 
bright XBF due to low metallicity. (The lower mass scale 
of Equation 7 could also be appropriate at solar metallic- 
ity if, in Section |4l we have overestimated the significance of 
stochasticity on XBF by an order of magnitude for any given 
SFR.) Such stellar masse s are consistent with the masses 
identified bv lSawala et al.l (|201ll ) for which feedback recipes 
in current galaxy models may be inadequate. More sophisti- 
cated investigation of this seems worthwhile. We stress that 
here we have estimated the mass scale at which significant 
XBF seems likely to be common in the greatest fraction of 
galaxies with that given mass. This is distinct from the con- 
tent of Section 14.11 in which we argue that a few percent of 
the lowest-mass halos should experience the greatest specific 
impact of XBF. 



6 DISCUSSION 

Feedback from XBs has properties that should not be ig- 
nored when modelling the formation and evolution of some 
galaxies. In particular, the long, slow energy input seems 
ideal to explain gas retention and episodic star formation 
in some low-mass dwarf galaxies. Even though metallicity is 
an important factor, the XBF probability distributions and 
the interaction of XBF and SNF are potentially even more 
important. The earliest population of galaxies, initially of 
very low metallicity and low absolute SFR, seem very likely 
to have been affected by ongoing XBF. 



6.1 The interplay between X-ray binary and 
supernova feedback 

Combining XBF and SNF suggests a rich variety of possi- 
bilities, especially given the broad probability distributions 
described in Section|4] We have explained that a few percent 
of halos with low absolute SFR might be expected to expe- 
rience a specific contribution from early XBF that is two or- 
ders of magnitude greater than the modal XBF input. Since 
the observed satellite galaxies may originate in only a sim- 
ilarly small percentage of haloes with low masses it seems 
reasonable to speculate on whether XBF could help some 
low-mass galaxies retain gas despite the effect of SNF. XBF 
should also apply to the metallicity-dependent disruption 
probabilities of globular clusters, which we intend to study 
in a future paper. XBF-driven expansion could at least de- 
crease the amount of gas within the virial radius at the point 
when SNe eject the bulk of the gas. 

Figure [5] illustrates three scenarios, all of which might 
occur in initially similar galaxies due to the stochastic varia- 
tion described in Section |4l Two of the scenarios occur when 
one of SNF or XBF is dominant; in either case the feedback 
might be strong enough to eject the gas reservoir from the 
galaxy. The intermediate case, when neither XBF or SNF are 
dominant, is potentially the most interesting. If XBF heat- 
ing is enough to stop star formation, but not to eject the gas, 
then the conditions into which SN explode are changed. The 
XBF-created hot bubble (s) seem(s) likely to aid the blow out 
of the SN ejecta (see, e.g.,[De Young & Hcckman 1994), and 
reduce the effective coupling of the SNF to the remainder of 
the gas reservoir. 

The balance between these three possibilities, and the 
quantitative conditions which define them, deserve further 
detailed work. For example, the development of chimneys 
and blow-out could be further aided if the distribution of 
the hot phase is strongly nonspherical, as might be expected 
for microquasars. Another intriguing detail which should be 
pursued is the difference in the ratio of energy and momen- 
tum involved in SNF and XBF. Moreover, in SS433 the dom- 
inant outflows of energy and mom entum are perpendicular 
to each other (JBlundell et al.ll200ll '). 



6.2 A scale length for supernova feedback? 

The effect that XBF pre-heating might have on SNF 
should be considered in all galaxies. Aware that pre-heating 
(by massive stars) might reduce the effectiveness of SNF, 
ICeverino fc KlvpirJ (|2009l ) suggested that those SNe pro- 
duced by runaway stars which explode in the colder gas out- 
side star-forming regions could dominate SNF. However, we 
note that this sets a natural scale-length. In some galaxies 
then stars may not be able to travel far enough to produce 
effe ctive SNF before leaving t h e gala xy. 

lEldridge. Langer fc ToutI (|201ll ) modelled populations 
of runaway supernovae and found that the most common 
class of supernova. Type IIP, ran away on average ~50 pc 
from their birthplace before exploding ( 48 pc at solar metal- 
licity, 54 pc at a metallicity of 0.004). If lCeverino fc KlvpinI 
(|2009f) are correct and the stars need to runaway to some 
extent in order to produce efficient SNF, then the average 
runaway distance might be significantly greater: consider- 
ing only the secondary stars which produce Type IIP su- 
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Figure 5. Schematic showing speculations on the potential interplay between XB and SN feedback. 



pernovae, i.e. selecting against those which don't run away, 
lEldridee. Langer &: Tout! ()201ll ) give a mean distance trav- 
elled of 190 pc at solar metallicity (and 280 pc at a metallic- 
ity of 0.004). The assumptions used by ICeverino fc KlvpinI 
(|2009f ) led to typical runaway distances between 10 and 100 
pc. (A runaway speed of 20 km s~^ for 10 Myr gives a dis- 
tance of «200 pc.) Runaway-star SNe might dominate SNF 
in more massive galaxies, yet decrease SNF in the lowest- 
mass dwarfs since more SNe take place outside the galaxy. 
A similar possibility emerges from a diflerent route: 
when the size of the hot bubble created by XBF dominates 
a galaxy, then again SNe seem less likely to be able to have 
the impact that they otherwise could. We note, as mentioned 
above, that the size of the X B-infiated nebula observed by 
IPakuU. Soria fc MotchI (|2010l 'l is « 300 pc, i.e. of compara- 
ble s ize to the lowest-mass Milk y Way satellite galaxies (see, 
e.g., iTolstov. Hill fc Tosi|[2009h . 



7 SUMMARY AND CONCLUSIONS 

Stellar physics should not be oversimplified if we desire accu- 
rate galaxy models; the physics of the baryonic component 
matters. We have explained why XB feedback could be im- 
portant, with characteristics perhaps significantly diflerent 
from SN feedback. XBF may well dominate some systems be- 
fore SNF acts, changing the properties of the galaxies into 
which SNF is later unleashed. XBF from young, luminous 
systems can still be important ~ 10* years after the for- 
mation of the relevant population, by which time the in- 
stantaneous star formation metallicity could have greatly 
increased. XBF from old populations might also be signifi- 
cant for keeping gas hot in some non-star-forming galaxies. 
The properties of XBF might stop star formation but 
still allow gas retention and multiple star-formation episodes 
in some very low-mass dwarf galaxies. The fraction of run- 
away massive stars might reduce the effectiveness of SNF for 
galaxies smaller than ~100 pc, even though those runaway 
SNe might increase the effect of SNF in massive galaxies. 



Furthermore, natural statistical variation in the balance be- 
tween energy input from XBs and SNe could help explain 
how most low-mass halos eject their gas but some satellite 
galaxies can retain gas during multiple SF episodes. The 
impact of XBF seems likely to be greatest in early-universe 
proto-galaxies, when low-mass halos are forming stars with 
low metallicity. In this case then episodic star formation 
might be expected in the early epochs of galeixy formation. 

Interplay between different forms of stellar feedback 
might increase the richness of galaxy evolution, especially 
for low-mass galaxies where the greatest statistical varia- 
tion in XBF is expected. The statistics of XBF in young 
populations also suggest that we should distinguish between 
the galaxy masses where the specific XBF should be largest 
in some individual cases (a tail of low-mass halos with low 
star formation rate and metallicity) and the galaxy masses 
where XBF seems most likely to impact a randomly-chosen 
halo (those which form ~ lO'^ M© of stars) . 

Probability distributions for metallicity-dependent 
XBF should be included in future galaxy formation mod- 
els, especially when studying dwarf galaxies and the earliest 
epoch of galaxy formation. 
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Taking the simplifying assumption of a uniform distri- 
bution of gas with mass Mg in a sphere of radius r, we can 
scale the the cooling luminosity Lcooi as approximately: 

A_22(^)7-^ ' (Al) 



where A_ 
10"^^ erg cm^ s 
ing. 



1036 erg s-~ ^ ^^"^^ [mZj VPC, 

cooling constant in units of 



is the 
^; this represents very efficient cool 



Al Keeping gas hot in old galaxies 

If we write the gas mass as a fraction of the stellar mass 
(Mt), such that Mg = gM*, then we can write the follow- 
ing for the ratio of the cooling to the XBF heating for old 
populations: 

^.10VA-..#^W" (A2) 

Llmxb Mq \pcy 

where Llmxb has been taken from Eq. [S] If we further 
use the mass-r adius relation derived for bright ellipticals in 
equation 12 ofJ Hasegan et al.l (|2005l ). we find: 

-0.845 

(A3) 



Lc, 



10 g A_22 



M, 



1012 M 



so a gas fraction of g ~ \/l0~6-3 ^ 0.1 per cent might be 
kept hot in such old massive ellipticals by the action of 
field LMXBs. In the very centre of a massive elliptical it 
seems reasonable to expect an aditional contribution from 
dynamically formed L MXBs (as in the inner bulge of M31, 
I Voss & Gilfanov|[2007h . as well as from LMXBs in globular 
clusters. Our cooling constant may also be too high for this 
hot gas regime, although our assumption of constant density 
will break down. 

We can also apply XBF to the other end of the ellipti- 
cal galaxy mass spectrum, i.e. to keeping gas hot in dwarf 
ellipticals. This is helped by the fact that they are over- 
sized compared to the mass-radius relation used above. For 
Mt = lO'' Mq and r = 1 kpc then g ^ 0.3 per cent. If the 
appropriate cooling constant is lower by a factor of 10 than 
our assumed value, then those those galaxies where the XBF 
luminosity distribution is enhanced by a factor of 10 above 
the mean (see Section [4]2]) could keep approximately 25 per 
cent of their baryonic mass in hot gas through the action of 
LMXBs. 



APPENDIX A: HEATING AND COOLING 
SCALING RELATIONS 

Here we provide approximate relations for the luminosity 
necessary to offset gas cooling, in order to illustrate the po- 
tential impact of XBF. Nonetheless, we wish to emphasize 
that kinetic outflows from some XBs suggest that the domi- 
nant effect of XBF might not be thermal. Furthermore, aver- 
aging XBF input over entire galaxies may also be misleading: 
the way that XBF locally restructures the ISM before SNF 
might be its most significant direct consequence (depending 
on the size of the galaxy in question) . 

In what follows we use the mean expected XBF energy 
inputs at solar metallicity, as presented in Section O i.e. 
Equations [l] and [3] (for the input to young and old popula- 
tions, respectively). 



A2 Star formation 

A similar cooling-to-heating ratio for XBF in young popu- 
lations at solar metallicity is: 

2 / \ -3 

(A4) 



icool _ 10-5 A / ^* 

« 5 X 10 A_22 --^ 



M, 



Moyr-i- 



We can simplify this by writing the star formation rate (M*) 
as Mgesp/rsF, where esF and tsf represent the efflciency and 
timescale with which gas is turned into stars. This gives: 

Ri5xl0 A_22^r^ — A5) 

LXBF eSF Mq \pcj 

i.e. a star formation timescale of 10 years and an efficiency 
of 5 per cent (or 10* years and 50 per cent) give a cooling-to- 
heating ratio of unity. Note that at low metallicity then not 
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only should XBF become stronger but cooling should also 
become less efficient. Furthermore, as discussed in Section 
m in many galaxies the XBF input varies from the numbers 
above simply via stochastic effects. 



